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\ The charged-current quasi-elastic scattering of muon neutrinos on oxygen target is com- 
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puted for neutrino energies between 200 MeV and 2.5 GeV using the relativistic distorted- 
^ - wave impulse approximation with relativistic optical potential, which was earlier successfully 

' applied to describe electron-nucleus data. We study both neutrino and electron processes 

' and show that the reduced exclusive cross sections for neutrino and electron scattering are 

<N : 

similar. The comparison with the relativistic Fermi gas model (RFGM), which is widely 
used in data analyses of neutrino experiments, shows that the RFGM fails completely when 
applied to exclusive cross section data and leads to overestimated values of inclusive and to- 
^ ' tal cross sections. We also found significant nuclear-model dependence of exclusive, inclusive 
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and total cross sections for about 1 GeV energy. 
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^ ■ PACS numbers: 25.30.-c 25.30. Bf, 25.30.Pt, 13.15.+g 



INTRODUCTION 
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' The field of neutrino oscillations has been rapidly developing from the observation of anomalies 

O ■ 

in cosmic rays [1] and solar [2] neutrino data to the cross checks of these anomalies [3, 4] and 
^ \ most recently to terrestrial confirmations of neutrino oscillation hypothesis (Kamland, K2K [5] and 

■ MINOS [6]). The next steps in this field would be the precision measurements of observed mass 

splitting and mixing angles and detailed experimental study of the neutrino mixing matrix. 

New, extremely intense neutrino beamlines are in operation or being planed. The data from 
these experiments will greatly increase statistics. In this situation, statistical uncertainties should 
be negligible compared to systematic uncertainties (ultimate precisions). An important source 
of systematic uncertainties is related to nuclear effects in neutrino interactions. Since nuclear 
targets are used as neutrino detectors, a reliable interpretation of neutrino data requires a detailed 
knowledge of energy and nuclear dependence of neutrino-nucleus {vA) cross sections. Apparently 
the uncertainties in neutrino cross sections and nuclear effects produce systematic uncertainties in 
the extraction of mixing parameters. 

Neutrino beams of high intensity cover the energy range from a few hundred MeV to several 
GeV. In this energy regime, the dominant contribution to neutrino-nucleus cross section comes 
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from quasi-elastic (QE) reactions and resonance production processes. Unfortunately, the cross 
section data in the relevant energy range are rather scarce and were taken on targets that are not 
used in neutrino oscillation experiments (i.e., water, iron, lead or plastic). 

A variety of Monte Carlo codes [7] developed to simulate neutrino detector response are based 
on a simple picture, referred to as Rclativistic Fermi Gas Model, in which the nucleus is described 
as a system of quasi-free nuclcons. Comparison with high-precision electron scattering data has 
shown that the accuracy of predictions of this model (inclusive cross sections) depends significantly 
on momentum transfer [8] . For inclusive nuclear scattering at sufficiently high momentum transfer 
(> 500 MeV/c) the RFGM describes general behavior of cross sections. However, the accuracy 
of a Fermi gas model becomes poor as momentum transfer decreases (see, e.g., [9]). Furthermore, 
this model does not account for the nuclear shell structure, and for this reason it fails when applied 
to exclusive cross sections. There arc other important effects beyond the RFGM: the final state 
interaction (FSI) between the outgoing nucleon and residual nucleus and the presence of strong 
short-range nucleon-nucleon (NN) correlations, leading to the appearance of high-momentum and 
high-energy components in the nucleon energy-momentum distribution in the target. In the calcu- 
lation of Ref. [10] within a plane-wave impulse approximation (PWIA) the nucleon-nucleon corre- 
lations were included using description of nuclear dynamics, based on nuclear many-body theory. 
It was shown that the Fermi gas model overestimates the total uA cross section by as much as 
20% at incoming neutrino energies of about 1 GeV. Neutral current and/or charged current (CC) 
neutrino-nucleus cross sections were studied within the relativistic distorted-wave impulse approx- 
imation (RDWIA) in Refs. [11-14] using a relativistic shell model approach. The implementation 
of the final-state interaction of the ejected nucleon has been done differently. A description of the 
FSI mechanisms through the inclusion of relativistic optical potential is presented in Refs. [11-13]. 
In Refs. [11, 12] important FSI effects arise from the use of relativistic optical potential within a 
relativistic Green's function approach. In Ref. [13], the final state interaction was included with 
and without the imaginary part of the optical potential (for inclusive cross section). A reduction of 
the total cross section of at least 14% was found at neutrino energies of 1 GeV. The relativistic op- 
tical potential and rclativistic multiple-scattering Glauber approximation were applied in Ref. [14] 
for the treatment of the FSI effects. Apart from relativistic and the FSI effects. Apart from rela- 
tivistic and FSI effects, other effects may be important in neutrino-nucleus reactions. In particular, 
Ref. [15-19] include long-range nuclear correlations (random-phase approximation) and FSI and 
Coulomb corrections in the calculation of i/^^C inclusive cross sections near threshold energy. 

In this paper, we compute the single-nucleon knockout contribution to the exclusive, inclusive. 
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and total cross sections of the charged-current QE (anti)neutrino scattering from ^^O using different 
approximations (PWIA and RDWIA) and the Fermi gas model. We employ the LEA code [20] 
developed for the calculation of contribution from Ip- and Is-state nucleons to cross sections 
in RDWIA. The LEA program, initially designed for computing of exclusive proton-nucleus and 
electron-nucleus scattering, was successfully tested against A(e,e'p) data [21-24], and we adopt 
this code for neutrino reactions. In the PWIA, the nuclear differential cross section are described 
in terms of a nuclear spectral function [25] , which includes contributions from nuclear shells as well 
as from the NN correlations. In our approach, the effect of the NN correlations in the oxygen 
ground state is evaluated in the PWIA using model nucleon high-momentum component [26, 27]. 
We propose a way to estimate the FSI effect on the inclusive cross sections in the presence of short- 
range NN correlations in the ground state. The aim of this work is twofold. First, we compute the 
RDWIA CC QE neutrino cross sections. Second, we test the RFGM against electron scattering 
data. 

The outline of this paper is the following. In Sec. II we present the formalism for the description 
of the charged-current lepton-nucleus scattering process. The RDWIA model is briefly introduced 
in Sec.III. Results of the numerical calculations are presented in Sec.IV. Our conclusions are sum- 
marized in Sec.V. In the appendix, we discuss the general Lorentz structure of the hadronic tensor 
and give expressions for the cross sections of neutrino exclusive scattering used in our analysis. 

II. FORMALISM OF QUASI-ELASTIC SCATTERING 

We consider electron and neutrino charged-current QE exclusive, 

l{ki) + A{pa) ^ l'{kf) + N{p^) + B{pb), (1) 

and inclusive, 

l{ki) + A{pA)^nkf) + X, (2) 

scattering off nuclei in a one-photon (W-boson) exchange approximation. Here I labels the incident 
lepton [electron or muon (anti) neutrino], and /' represents the scattered lepton (electron or muon). 
Figure 1 defines our conventions for the kinematical variables, where ki = (sj, fcj) and kf = (ej, kf) 
are initial and final lepton momenta, pA = {£a,Pa)j ^^id pB = {£b-,Pb) &re the initial and final 
target momenta, Px = {ex,Px) is ejectile nucleon momentum, q = {LO,q) is the momentum transfer 
carried by the virtual photon (W-boson), and = —q^ = — uP' is the photon (W-boson) 



reaction plane 



FIG. 1: (Color online) Kinematics for the quasi-elastic lepton-nucleus scattering process, 
virtuality. Normalization of states is given by 

NM\P'^='^^SHP^-P'i), 

where Ni = ra/e for massive particles, or Ni = l/2e for massless leptons. 

A. Quasi-elastic lepton-nucleus cross sections 

In the laboratory frame, the differential cross section for exclusive electron (c*^') and 
(anti) neutrino {<j'^'^) CC scattering can be written as 



(ffjCC \Px\^x l^/l G 

dsfdrifdexdnx ~ (27r)5 ~i 2 ^f"" 



where is the solid angle for the lepton momentum, 0,^ is the solid angle for the ejectile nucleon 
momentum, a ~ 1/137 is the fine-structure constant, G ~ 1.16639 xlO^^^ MeV^^ is the Fermi 
constant, 6c is the Cabbibo angle (cos 6'c ~ 0.9749), L^'^ is the lepton tensor, W^Ji^ and wj^i^^ are 
correspondingly the electromagnetic and weak CC nuclear tensors which will be discussed below. 

For exclusive reactions in which only a single discrete state or narrow resonance of the target is 
excited, it is possible to integrate over the peak in missing energy and obtain a fivefold differential 
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cross section of the form 

I, _ ^,2 



cc 



p JPjgx \kf\G^ COS^ gc ^(cc) ^mKcc) ^ 



defdftfdn^ (27r)5 2 '''' 

where i? is a recoil factor 



i? = y dexS{ex +£b -(^ - "t-a) 



Px ■ Px 



(4b) 



(5) 



£x is the solution to the equation Sx + £b — fnA — a; = 0, where es = y + p^, = q — p^ 
and m^i and are masses of the target and recoil nucleus, respectively. Note that the missing 
momentum is = Px ~ 

The lepton tensor can be written as the sum of the symmetric Lg^ and antisymmetric L'^ 
tensors 

L'^'^ = + (6a) 

L^'^ = 2 (A;f fc^ + fc^^/fc^ - g^^kikf^ (6b) 

l:^'^ = h2iei''"'P(ki)^{kf)p, (6c) 

where /i is +1 for positive lepton helicity and —1 for negative lepton helicity, and e^^°^^ is the 
antisymmetric tensor with e'^^^^ = — eoi23 = 1- For the scattering of unpolarized incident electrons, 
/,Mi'(eO Q^\j jjg^g ^]jg symmetric part (6b) and the (anti) neutrino tensor L^^^'^'^^ involves both the 
symmetric and the antisymmetric parts. Assuming the reference frame, in which the z axis is 
parallel to the momentum transfer q = ki — kf and the y axis is parallel to fej x fej, the symmetric 
components iPg, L^^, L^g and the antisymmetric ones L^, L^^, L^^, as well as those obtained from 
them by exchanging their indices, vanish. The electromagnetic and the weak CC hadronic tensors, 

ieX) [cc] 

W'^u' and W^r, are given by bilinear products of the transition matrix elements of the nuclear 
electromagnetic or CC operator jj^''^^'^'^^ between the initial nucleus state \A) and the final state 
\Bf) as 

Wg)^'^^) = Ji^')('=^)|^)(^|4'^')('-)t|s^,^^)^(,^ +u-ex- SB,), (7) 

/ 

where the sum is taken over undetected states. 

In the inclusive reactions (2) only the outgoing lepton is detected, and the differential cross 
sections can be written as 



defdflf £i g4 t^'^ 
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defdVLf (27r)2 ei 

where W^'^ is inclusive hadronic tensor. A general covariant form of the hadronic tensors and the 
results of their contractions with the lepton tensors are given in Appendix A for exclusive lepton 
scattering (1). Combining Eq.(4a) with Eq.(A3) and Eq.(4b) with Eq.(A7) we obtain the exclusive 
lepton scattering cross sections in terms of response functions 



defdnfdn " ^^^^^(^^^i + + ^ltR^l? cos + VttR!^} cos 24>) , (9a) 



G cos OcSf\kf\R[vQRo + VtRt + VTTRTTCOS2(j) + VzzRz 



dsfdQfdn^ (27r)5 

+ {vxzRxz - vqxRqx) cos (p - vqzRqz + h [vyz{Ryz sin (j) + Ryz cos (j)) 

- voy (R^y sin 4> + Roy cos 0) - VxyRxy] }, (9b) 



where 



= .-2.^,4;,. (10) 



cos^ e/2 
iefsh?^ 

is the Mott cross section. The response functions Ri depend on the variables Q^,uj, \Px\, and 9x- 
Similarly, the inclusive lepton scattering cross sections reduce to 

^^^^^^ = (27r)2 "" ^/l^/K^o-Ro + vtRt + VzzRzz - vozRoz - hvxyRxy), (lib) 

where the response functions now depend only on and uj. 
It is also useful to define a reduced cross section 

'"'^^'^ = dsf^fdnx^^'"'''' ^^^^ 

where = Rpx6x/{27r)^ and K"^^^ = Rp^Sx/ {2tt)^ are phase-space factors for the electron and 
neutrino scattering, the recoil factor R is given by Eq.(5), and ctjjv is the corresponding elementary 
cross section for the lepton scattering from the moving free nucleon. 



B. Nuclear current 



Obviously, the determination of the response tensor requires the knowledge of the nuclear 
current matrix elements in Eq.(7). We describe the lepton-nucleon scattering in the impulse ap- 
proximation (lA), assuming that the incoming lepton interacts with only one nucleon, which is 
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subsequently emitted. The nuclear current is written as the sum of single-nucleon currents. Then, 
the nuclear matrix element in Eq.(7) takes the form 

{p,B\J^'\A) = exp(it-r)*^"\p,r)r'^$(r), (13) 

where V^^ is the vertex function, t = ebq/W is the recoil-corrected momentum transfer, W = 
{jriA + fjS)'^ — cp" is the invariant mass, $ and ^'^"^ are relativistic bound-state and outgoing 
wave functions. 

For electron scattering, most calculations use the CC2 electromagnetic vertex function for a free 
nucleon [32] 

= Ff\Q')^^ + ia^-^F^^\Q\ (14) 

where a^"" = i['y^' Fy^^ and fJ^'^ are the Dirac and Pauli nucleon form factors. Because the 
bound nucleons are off shell, the vertex in Eq.(13) should be taken for the off-shell region. We 
employ the de Forest prescription for off-shell vertex [32] 

where q = (e^- — E, q) and the nucleon energy E = ^ m? + (p^ — qY is placed on shell. We use 
the approximation of [33] on the nucleon form factors. The Coulomb gauge is assumed for the 
single-nucleon current. 

The single-nucleon charged current has the V—A structure J'*('^) = + J^. For a free nucleon 
vertex function Vi^^'^^ = F^ -|- F^, we use the CC2 vector current vertex function 

F^ = FviQ^h'' + ia^''^FM{Q% (16) 

and the axial current vertex function 

F^^ = FA{Q^)r^5 + Fp(Q')g'^75. (17) 

The weak vector form factors Fy and Fm are related to the corresponding electromagnetic ones 
for proton Fj^^^'' and neutron F^^^^ by the hypothesis of the conserved vector current (CVC) 

F^ = F^f-FS^■ (18) 
The axial Fa and psevdoscalar Fp form factors in the dipole approximation are parameterized as 

where -Fa(O) = 1.267, m,r is the pion mass, and Ma — 1.032 GeV is the axial mass. We use the de 
Forest prescription for off-shell extrapolation of F^^^'^). Similar to the electromagnetic current, the 
Coulomb gauge is applied for the vector current Jy- 
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III. MODEL 

In Ref. [34], a formalism was developed for the A{e,e'N)B reaction that describes channel 
coupling in the FSI of the N + B system. According to Ref. [34] , a projection operator P for 
model space was introduced. In the independ particle shell model (IPSM), the model space for 
^^0(e, e'iV) consists of lsi/2, lP3/2; ^^'^ lPi/2 nucleon-hole states in -"^^N and -^^O nuclei, for a total 
of six states. The lsi/2 state is regarded as a discrete state even though its spreading width is 
actually appreciable. For single nucleon knockout, the parentage expansion of the target ground- 
state can be written as 

P*o = 5Ic/3t'/'/3t*7' (20) 

where cpj is a parentage coefficient and (pp^ is an overlap wave function for removal of a nucleon 
with single-particle quantum number /? while leaving the residual nucleus in the state Assuming 
that the overlap wave functions are described by the Dirac equation, they can be represented by a 
Dirac spinor of the form 

^/J7 = I I • (21) 

Similarly, for the scattering state 

P*W=^VS^$/3 (22) 

is an incoming wave function of the N + B system containing an incident plane wave in the 
channel a and outgoing spherical waves in all open channels /? for B{N, N')B' reaction. The Dirac 
representation of distorted spinor wave functions is 

^iV = r , (23) 

where 




is the asymptotic wave function for channel a normalized to unit flux, and = y/k^+'m^ is the 
channel energy in the barycentric frame (the rest frame of residual nucleus B). 

Working in coordinate space, we can write the matrix elements of the current operator (16) for 
single-nucleon knockout leaving the residual nucleus in asymptotic channel a as follows 

{j),Ba\J^'\A) = cp^ J d^rexp{it-r)0J\rmb) 
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x(rm6|f'*|rm;)(rm|,|<^;3^). (25) 
Matrix elements of the single-nucleon current can be expressed in the block-matrix form 

f'' = , (26) 




where each of the elements {rmb\r'^y k^i^) is a 2 x 2 spin matrix, while A = {+, — } and A' = {+, — } 
are for the upper (+) and lower (— ) Dirac components. Let 



(27) 



be the bound state overlap wave function and 



(Viy|rm,) = A^J 1 (28) 

be the Dirac adjoint of time-reversed distorted waves. 

For the sake of application to cross section calculations, we consider the relativistic bound-state 
functions within the Hartree-Bogolioubov approximation in the a-co model [35]. In the mean-field 
approximation, the meson field operators are replaced by their expectation values. The upper and 
lower radial wave functions in the partial-wave expansion for bound-state wave functions satisfy 
the usual coupled differential equations 

(J^ + '^^ypy{r) = [E^ + m + S^ir)-V^{r)]G0^ir), (29a) 

{'i ~ ^^-'^''^ =[-E, + m + S,{r) + V^{r)] Fp^{r), (29b) 

where Sj and Vy are spherical scalar and vector potentials, and j-y = |K-y| — 1/2 is the total angular 
momentum. Note that these potentials generally depend on the state of the residual nucleus that 
is marked by subscript 7. The radial wave functions are normalized as 

J drr^[\Fp^\^ + \Gp^\^) =1. (30) 

The missing momentum distribution is determined by the wave functions in momentum space 

Fp'riP) = J dr r^ji^{pr)Fp^ir), (31a) 
,{p) = J dr r^ji>{pr)G^^{r), (31b) 



G 
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where ji{x) is the Bessel function of order I and l'^ = 2jj — Ij. If only a single state of residual 
nucleus is considered, or if relativistic potentials S and V weakly depend on the state 7 of residual 
nucleus, the relativistic momentum distribution can be written in terms of Eq.(30) as 



(32) 



In this work, the current operator CC2 and the bound-nucleon wave functions [36] (usu- 
ally referred to as NLSH) are used in the numerical analysis. Note that the calculation of the 
bound-nucleon wave function for state includes the incoherent contribution of the unresolved 
251/2^^5/2 doublet. The wave functions for these states were taken from the parameterization of 
Ref. [37]. We use also the following values of normalization factors Sa = \ca\^ relative to the full 
occupancy of ^^O: ^(Ipg/a) = 0.66, ^(Ipi/a) = 0.7 [23], and S{lsy2) = 1- 

The distorted wave functions are evaluated using a relativized Schrodinger equation for upper 
components of Dirac wave functions. For simplicity, we consider a single-channel Dirac equation 



[a-p + Pirn + S)]^; = {E - V)^;, 



(33) 



where 



^(r) 




(34) 



is the four-component Dirac spinor. Using the direct Pauli reduction method [38, 39], the system 
of two coupled first-order radial Dirac equations can be reduced to a single second-order equation 



(35) 



where ^ is a two-component Pauli spinor. Here k is the relativistic wave number, /x is the reduced 
mass of the scattering state, and 



1^ 



E 



+ u 



D 



2/1 



U 



LS 



1 d 
2r^D dr 
1 D' 



2/xr D ' 



(36a) 
(36b) 
(36c) 
(36d) 



E + m 

where D' = dD/dr, and D{r) is known as the Darwin nonlocality factor, and U'-^ and U^^ are the 
central and spin-orbit potentials. The upper and lower components of the Dirac wave functions 
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axe then obtained using 



i,+ = D^l^i, {37a) 

= g + J+g-F 

Assuming a similar relationship for the coupled-channel case, i.e., 

the lower components of the radial wave functions in the partial-wave expansion for distorted waves 
(31) can be approximated as 



(39) 



We use the LEA program [20] for the numerical calculation of the distorted wave functions with 
the EDADl SV relativistic optical potential [40]. This code employs an iteration algorithm to solve 
the relativized Schrodinger equation. 

A complex relativistic optical potential with a nonzero imaginary part generally produces an 
absorption of flux. For the exclusive channel, this reflects the coupling between different open 
reaction channels. However, for the inclusive reaction the total flux must conserve. Currently 
there is no fully consistent solution to this problem, and different approaches are used. The Green's 
function approach, where the FSI effect in inclusive reactions is treated by means of a complex 
optical potential and the total flux is conserved, is presented in Refs.[ll, 41]. To demonstrate 
the effect of the optical potential on the inclusive reactions, the results obtained in this approach 
were compared with those obtained with the same potential but with the imaginary part set to 0. 
It was shown that the inclusive CC neutrino cross sections calculated with only the real part of 
optical potential are almost identical to those of the Green's function approach [11, 12]. A similar 
approximation was used also in Ref. [13] to study the FSI effect on the inclusive cross section. In 
this work, in order to calculate the inclusive and total cross sections, we use the approach in which 
only the real part of the optical potential EDADl is included. Then the contribution of the Ip and 
Is states to the inclusive cross section can be obtained by integrating the exclusive cross sections 
(11) over the azimuthal angle ^ and missing momentum, that is, p„ 



apm — i— r-^c X 



= d<f> 

Jo J Pmin 



d£fd^f J RBWiA Jo Jpmin Px\q\ \d£ fd^fd^x J rbwia 
where pm = Px = \Px\^ Pm = Px - g> and 

cos 9x = + ~ , (40a) 
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Rc = l + ^{pl + q'-pl)- (40b) 

The integration hmits Pmin and Pmax are given in Ref. [27] . The effect of the FSI on the inclusive 
cross section can be evaluated using the ratio 

= (5$i7)pwA' 

where {d^a/dsfd^f) p.^^^ is the result obtained in the PWIA. 

According to data from the Thomas Jefferson National Accelerator Facility (JLab) [23], the 
occupancy of the IPSM orbitals of -"^^O is approximately 75% on average. In this paper, we assume 
that the missing strength can be attributed to the short-range NN correlations in the ground state. 
To estimate this effect in the inclusive cross sections, we consider a phenomenological model. This 
model incorporates both the single particle nature of the nucleon spectrum at low energy and high- 
energy and high-momentum components due to NN correlations. The high-momentum part Phm 
of the spectral function is determined by excited states with one or more nucleons in continuum. 
The detailed description of this model is given in Refs.[26, 27]. 

In our calculations the spectral function Phm incorporates 25% of the total normalization of 
the spectral function. The FSI effect for the high-momentum component is estimated by scaling 
the PWIA result {d^a/d£fdQ,f)HM with A{ef,Qf) function (41). Then the total inclusive cross 
section can be written as 

d^a ( d?a \ , , f cP'cr \ 

^(^/'^/)b73^ • (42) 



dsfd^f \d£fd^fJ^Yi^ip^ \d£fd^fJmA 
More details about calculation of the {d^a/d£fdQ,f)uM can be found in Ref. [8]. 



IV. RESULTS 



The LEA code was successfully tested against A{e, e'p) data [22-24]. In Ref. [23] the uncertainty 
in the normalization factors Sa was estimated to be about ±15%. For illustration. Fig. 2 shows 
the measured JLab [23] and Saclay [42] differential cross sections for the removal of protons from 
the Ip shell of ^^O as functions of missing momentum pm as compared with LEA code calculations. 
The reduced cross sections together with Saclay [43] and NIKHEF [44] data are shown in Fig. 3. 
It should be noted that negative values of Pm correspond to <j) = it and positive ones to (j)=Q. The 
cross sections were calculated using the kinematic conditions with the normalization factors of data 
examined [23]. Also shown in Figs. 2 and 3 are the results obtained in the PWIA and RFGM 
(with the Fermi momentum pF=22h MeV/c, binding energy e=27 MeV and including the Pauli 



13 




FIG. 2: (Color online) Calculations compared with measured differential exclusive cross section data for the 
removal of protons from the Ip shell of ^^O as a function of missing momentum. Upper panels: JLab data 
[23] for beam energy £;bcam=2.442 GeV, proton kinetic energy Tp^427 MeV, and Q'^=0.8 GeV^. Lower 
panels: Saclay data [42] for £;bcam=580 MeV, Tp=160 MeV, and Q'^=0.3 GeV^ 

blocking factor). Apparently the PWIA and RFGM overestimate the values of the cross sections, 
because the FSI effects are neglected. Moreover, the RFGM predictions are completely off of the 
exclusive data. This is because of the uniform momentum distribution of the Fermi gas model. 
The reduced cross sections for the removal of nucleons from Ip shell in ^^0(e, e'p)^^N, 
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FIG. 3: (Color online) Calculations compared with measured reduced exclusive cross section data for the 
removal of protons from the Ip shell of ^^O as a function of missing momentum. Upper panels: Saclay 
data [43] for beam energy £'bcam=500 MeV, proton kinetic energy Tp— 100 MeV, and Q^=0.3 GeV^. Lower 
panels: NIKHEF data [44] for £'bcam=521 MeV, Tp=96 MeV, and is varied. 

^^0{i^, H~p)^^0, and ^^0(P, /i^n)^^N reactions are shown in Fig. 4 as functions of pm together 
with Saclay [43] and NIKHEF data. There is an overall good agreement between calculated cross 
sections, but the value of electron cross sections at the maximum is systematically higher (less than 
10%) than (anti)neutrino ones with the exception of the lpi/2 state for Saclay kinematics. The 
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FIG. 4: (Color online) Comparison of the RDWIA electron, neutrino and antineutrino reduced cross sections 
for the removal of nucleons from the Ip shell of ^^O for Saclay [43] and NIKHEF [44] kinematic as functions 

of Pm- 

small difference between neutrino and antineutrino reduced cross sections is due to the difference 
in the FSI of proton and neutron with the residual nucleus. 

The differential and reduced electron and (anti)neutrino exclusive cross sections for the removal 
of nucleons from Ip and Is states were calculated for JLab and Saclay [42] kinematics. The results 
are shown in Fig. 5 together with the RFGM calculations. There is a good agreement between all 
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FIG. 5: (Color online) Comparison of the RDWIA and the RFGM calculations for electron, neutrino and 
antineutrino reduced (left panels) and differential (right panels) cross sections for the removal of nucleons 
from Ip and Is shells of ^^O as functions of missing momentum. The cross sections were calculated for the 
JLab [23] and Saclay [42] kinematics. In the left panels, the RDWIA calculations are shown for electron 
scattering (dashed-dottcd line) and neutrino (dashed line) and antineutrino (dotted line) scattering; and 
the RFGM results are shown for the reduced cross sections (solid line). In the right panels, the RFGM 
calculations are shown for the neutrino (solid line) and antineutrino (dashed line) differential cross sections; 
and the RDWIA results are shown for the neutrino (dashed-dotted line) and antineutrino (dotted line) 
differential cross sections. 
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cross sections calculated in the RDWIA for JLab kinematics. The difference between the electron 
and (anti)neutrino reduced cross sections calculated for Saclay kinematics is less than 10%. This 
can be attributed to Coulomb distortion upon the electron wave function which is usually described 
in the effective momentum approximation (EMA) [45]. In the EMA, the electron Coulomb wave 
function is replaced by a plane wave with effective momentum whose value is larger than the 
value of electron momentum at infinity, because of Coulomb attraction. This effect weakens as 
the beam energy increases, and for this reason this effect is more significant at Saclay kinematics 
(^^beam = 500 MeV) than at JLab kinematics (£^beam = 2442 MeV). Note that the RFGM results 
demonstrate absolutely different behavior. 

To test our approach, we calculated the inclusive ^^0(e, e') cross sections and compared them 
with SLAC data [46] and Prascati data [47] . Figures 6 and 7 show measured inclusive cross sections 
as functions of energy transfer, or the invariant mass W as compared with the RDWIA, PWIA, 
and RFGM calculations. We note that relative to the PWIA results, the generic effect of the FSI 
with the real part of the optical potential is to reduce the cross section value around the peak and 
to shift the peak toward the lower value energy transfer. The inclusion of the high-momentum 
component increases the inclusive cross section in the high-energy transfer region and improves 
the agreement with data. For the RDWIA results, the difference between the calculated and 
measured cross sections at the maximum are less than ±10%, with the exception of Prascati data 
for Eg = 700 MeV. For the RFGM results, these differences decrease with |g| from about 22% at 
|q| ~ 330 MeV/c down to ~ 2% at \q\ ~ 640 MeV/c. These results demonstrate a strong nuclear- 
model dependence of the inclusive cross sections at low momentum transfer. This dependence 
weakens as |q| increases, almost disappearing at |q| > 500 MeV/c. The results for {e,e'N) channel 
indicate that at least 50% of the inclusive cross section can be attributed to the single-step nucleon 
knockout. 

The inclusive neutrino and antineutrino cross sections for energies = 300, 500, 700, and 1000 
MeV are presented in Figs. 8 and 9, which show da/dEfj, as a function of muon energy. Here, the 
results obtained in the RDWIA with the real optical potential (RDWIA ROP) are compared with 
the inclusive cross sections calculated in the PWIA, RFGM, and RDWIA with complex optical 
potential (RDWIA EX). The cross section values obtained in the RFGM are higher than the ones 
obtained within the RDWIA ROP. For neutrino (antineutrino) cross sections in the region close to 
the maximum, this discrepancy is about 35%(60%) for E^ = 300 MeV and 30%(40%) for E^ = 1000 
MeV. The contribution of {v^ixN) channels to the inclusive cross sections is about 60%. 

The total cross sections a {En) together with data [48, 49] are presented in Fig. 10 as functions 
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FIG. 6: (Color online) Inclusive cross section vs energy transfer uj or invariant mass W for electron scattering 
on 1*^0. The data are from Ref.[46] (SLAC, filled circles) and Ref.[47] (Frascati, filled triangles). SLAC data 
are for electron beam energy Ee—54Q, 730 MeV and scattering angle 0e=37.1°. Frascati data are for Ee—540 
MeV and 6'e=37.1°, Ee—700, 880 MeV and 9e—32°. As shown in the key, cross sections were calculated with 
the RDWIA, PWIA, RFGM and RDWIA with complex optical potential (EX). 

of the incident neutrino energy. The upper panel shows the total cross sections for ^^O {u, 
reaction calculated in the RDWIA with the real part of EDADl potential, and the lower panel 
shows the total cross sections for the ^^0(f^, /x~p) channel. Also shown are the results obtained in 
Refs. [11, 13] with the NLSH bound nucleon wave functions, dipole approximation of the nucleon 
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FIG. 7: (Color online) Same as Fig. 6, but the data are from Ref.[47] for electron beam energy i?e=1080, 
1200, and 1500 MeV and scattering angle 6'e=32°. 

form factors, EDADl optical potential and neglecting the NN correlation contributions. The cross 
sections are scaled with the number of neutrons in the target. 

Our values of (z^, fJ-^p)] cross sections are systematically larger than those from Ref. [11]. 

The discrepancy increases with energy from about 17%(7%) for Ei, = 300 MeV up to 28%(20%) 
for = 1000 MeV. On the other hand, our cross sections are lower than those from Ref. [13], and 
the discrepancy decreases with energy from 37%(15%) for Ei, = 300 MeV upto 15% (7%) down 
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FIG. 8: (Color online) Inclusive cross section vs the muon energy for neutrino scattering on -'^^O and for the 
four values of incoming neutrino energy: 0.5, 0.7, and 1 GeV. 

to El, = 1000 MeV. To study the NN correlation effect, we calculated the total cross sections 
without the high-momentum contribution, i.e., with Sa = 1 for all bound nucleon states, similar 
to Refs [11, 13]. The results are shown in Fig. 10. Apparently, the A^A^ correlation effect reduces 
the total cross section. The difference between the results obtained with and without the high- 
momentum component contribution decreases with neutrino energy from about 20% for = 200 
MeV down to ~ 8% for Ei, = 1000 MeV. Moreover, in this case the agreement with the result of 
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FIG. 9: (Color online) Same as Fig. 8, but for antineutrino scattering. 

Ref.[13] is good, and the discrepancy is less then ±6% for E,^ > 300 MeV. 

The neutrino and antineutrino total cross sections calculated up to neutrino energy 2.5 GeV 
are shown in Fig. 11 together with data of Refs. [48-51]. Also shown are the results obtained in the 
RFGM and PWIA as well as the contribution of the exclusive channels to the total cross sections. 
The cross sections are scaled with the neutron/proton number in the target. The ratio between the 
neutrino cross sections calculated in the RFGM and RDWIA ROP decreases with neutrino energy 
from about 1.5 for = 300 MeV to w 1.18 for E^ = I GeV and down to w 1.05 for E^ = 2.4 
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FIG. 10: (Color online) Total cross section for the CC QE scattering of muon neutrino on -'^^O as a function 
of the incoming neutrino energy. The RDWIA results with the real part of optical potential (upper panel) 
and complex optical potential (lower panel) are shown together with calculations from Meucci et al. [11] and 
Maieron et al. [13]. The results obtained in this work were calculated with and without the contribution of 
the high- momentum component. For comparison, data for the D2 target are shown from Refs.[48, 49]. 

GeV. For the antineutrino cross sections, this ratio is about 2.7 for Ei, = 300 MeV, 1.3 for = \ 
GeV, and 1.1 for E^ = 2.4 GeV. 

It follows from the comparison of the PWIA and RDWIA results that the FSI effects reduce 
the total cross section. For the neutrino interactions, this reduction is about 16% for Ey = 300 
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FIG. 11: (Color online) Total cross section for CC QE scattering of muon neutrino (upper panel) and 
antineutrino (lower panel) on ^^O as a function of incoming (anti) neutrino energy. Data points for different 
targets are from Refs. [48-51]. 

MeV and decreases slowly to 10% for Ei, = 2.4 GeV. The reduction of the antineutrino cross 
section is about 38% for = 300 MeV and « 15% for E^, = 2.4 GeV. We, therefore, observe the 
weakening of FSI effect in total cross sections with the increase of energy transfer, in accordance 
with the calculation of Ref.[52]. The contribution of the exclusive channels is about 60%. The 
results presented in Fig. 11 show significant nuclear-model dependence for energy less than 1 GeV. 
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V. CONCLUSIONS 

In this paper, we study electron and CC quasi-elastic (anti)neutrino scattering on the oxygen 
target in different approximations (PWIA, RDWIA, RFGM) placing particular emphasis on the 
nuclear-model dependence of the results. In RDWIA, the LEA program, adapted to neutrino 
interactions, was used to calculate the differential and reduced exclusive cross sections. This 
approach was earlier applied to electron-nucleus scattering and successfully tested against data. 
We found that the reduced cross sections for (anti)neutrino scattering arc similar to those of 
electron scattering, and the latter are in a good agreement with electron data. In calculating the 
inclusive and total cross sections, the imaginary part of relativistic optical potential was neglected 
and the effect of NN correlations in the target ground state was taken into account. This approach 
was tested against electron-oxygen inclusive scattering data; there was overall agreement with the 
data, with the differences between calculated and measured cross sections in the peak region less 
than 10%. For neutrino interactions the FSI effect reduces the total cross section by about 30% 
for E^=200 MeV compared to PWIA and decreases with neutrino energy down to 10% at 1 GeV. 
The effect of NN correlations reduces the total cross section by about 15% at £'i/=200 MeV and 
also decreases with neutrino energy down to about 8% at 1 GeV. 

We tested the RFGM against electron-oxygen scattering data and found that this model does 
not reproduce the exclusive cross section data. The RFGM also leads to an overestimated value 
of the inclusive ^^0(e, e') cross section at low momentum transfer. The discrepancy is about 20% 
and decreases as momentum transfer increases. The values of the (anti)neutrino cross sections 
calculated in this model are also higher than the corresponding values in the RDWIA approach. 

We conclude that the data favor the RDWIA results. This indicates that the use of RDWIA in 
Monte Carlo simulations of neutrino detector response would allow one to reduce the systematic 
uncertainty in neutrino oscillation parameters. 

Acknowledgments 

The authors greatly acknowledge communications with J. J. Kelly whose LEA code for nucleon 
knockout by electron scattering was adapted in this work for neutrino interactions. This work 
was partially supported by Russian Foundation for Basic Research, Project Nos. 06-02-16353, 
06-02-16659, and 05-02-17196. 



25 



APPENDIX A: HADRONIC TENSOR AND CROSS SECTION OF EXCLUSIVE 
ELECTRON AND NEUTRINO SCATTERING 

A general structure of the hadronic tensor can be derived from the requirements of Lorentz 
invar iance, parity, and time reversal symmetries. For unpolarized nucleon and a nucleus in the 
final state, this tensor must be constructed from three linearly independed four- vectors q, Px, 
and PA, the scalars that can be constructed from them, the second-rank metric tensor g^t^,, and 
completely antisymmetric tensor e^^Q,^. Generally, because of the final state interaction effects, 
the scattered flux at infinity involves complicated asymptotic configurations, and the time reversal 
symmetry does not constraint the form of the nuclear tensor for the exclusive reactions [28-30] . 



1. Electron scattering 

For electron scattering, the electromagnetic current conservation requires q^W^^ = W^^q^ = 0. 
Taking into account the parity conservation, the nuclear tensor can then be written as a sum of 
symmetric, Wg^^^^\ and antisymmetric, W^^^^^\ parts [31] 



^/.^eO ^ ^^Hel) + ^M-(eO ^ ^ ^^^^ 

W^'^iel) ^ ^(eO ^ W^^^VxP^ + W^'VaPA + Wt'^ mi + &a\ (Alb) 

WT^''^ = W^^'\p^xP'a-pIp'a)^ (Ale) 



where 



r"=5''" + ^, (A2a) 
P^=P^ + ^9^ (A2b) 
^a=p'a + '^<i'- (A2c) 

In target rest frame, the coordinate system is chosen such that the z axis is parallel to the momen- 
tum transfer q = kj — kj and the y axis is parallel to kj x kj, and the components of the four- 
vectors are kf = {sf, |k/| sin6cos<^, |k/| sin^sini^, |k/|cos6), q = (u;,0,0, |q|), pA = (m^, 0,0,0), 
Px = (^s) jPij sin^jrcos (?!), |pj;| sin^j; sin^, jp-rj cos^j;)) where 6, tp are lepton scattering angles and 
Ox, (j) the outgoing nucleon angles. 

The lepton tensor for unpolarized electron scattering is symmetric, and therefore the result of 
contraction of the electron and nuclear response tensors reduces to the form 
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Lg) W^"^^^^ = AsiE f cos^ - {VlR'i^^ + VtR^t^ + VltR^t cos <P + VttRtt cos 2(^) , (A3) 

where 

Vl = QVq\ (A4a) 

^T = ^ + tan2^, (A4b) 



2g2 



Ftt = (A4d) 



are the electron couphng coefficients, and 

^(eO^^OO(eO^ (A5a) 

COS = - + ^ (A5c) 

COS 2(t> = - W^s'^^^^) , (A5d) 

are four independ response functions, which describe the electromagnetic properties of the hadronic 
system. 

2. Neutrino scattering 

In weak interactions, the weak current and parity are not conserved. Therefore, a general 
nuclear tensor can be written as 



^M^(cc)^ W^^' + W'^^, (A6a) 

+ W^Ip^Xq^'+pW) + W,{jr,p\+plp^Xl (A6b) 
WT = Wsip^^q" - Plq^) + - PaQ^ + W^P^^p^X - pIp\) 

+ Wxxe^'^-PqrP^p + Wx2e^''^PqrPAp + W^ze^^^^fp^rPAp- (A6c) 

Note that because of Hermicity of W^^'^"^ , each term of W*^ must be real, while each term of W*^ 
must be imaginary, and U^JW^^^"^ is real. The result of contraction of the lepton and nuclear 
tensors can be written as 
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where 



^(cc)^mKcc) = L^^^w^- + i^^^W^r = 2£i£/{i;oi?o + vtRt + vttRtt cos 20 + V,,R,, 
+ {vxzRxz - vqxRqx) cos (j) - vqzRqz + h[vyz{Ryz sin (f) + Ryz COS 

- voyiR'oy sin 'Z' + -Roy cos 0) - 'Uiyi^ij/]}, (A7) 

= 1 + /3cos6', (A8a) 

.^ = l-/3cos^ + £^^^^^, (A8b) 

ei/?|fe/-| sin^6' , 

VTT = -^^^^ , (A8c) 

2 

^^0^ = A(l + /3cose) + ^, (A8d) 

i,^^ = 1 + /? cos ^ - 2^^^^ sin^ e, (A8e) 

, , B sin . , „ „s 

■^sm9[{ei + £f)u} + m'{)], (A8g) 
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2 



{1- Pcos9)- (A8h) 

I 

t;2;2 = /3|^sin^, (A8i) 

voy = Psm9, P=\kf\/ef, (A8j) 



are neutrino coupling coefficients, and 



Ro = Wl°, (A9a) 

= W^"" + (A9b) 

i?Tr cos 20 = ^ - W^y, (A9c) 

i?o. = W^^^ + Wf, (A9d) 

i?^, = Wf, (A9e) 

i?o^cos0 = + (A9f) 

i?^,cos0 = W^|^ + PF|^ (A9g) 

Rxy = i{W^'-Wf), (A9fl) 

i?^^ sin + i?^/^ cos = z ( - W^^) , (A9i) 
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E^Qy sin cf> + Roy cos (t> = i - Wf^ , (A9j ) 

are ten independ response functions which describe the weak properties of the hadronic system. 

In the absence of FSI effect (plane-wave hmit) the nucleon flux conserves in exclusive reaction. 
For this reason, the time reversal symmetry of operators and states provides an additional con- 
straint on the Lorenz form of the antisymmetric part of nuclear tensor (A6c), in particular, the 
structures like a'^b'^ — a^h^ (with a and h four-momenta) vanish. Then we have 

Wr = W^^e^'''^PqrP.p + Wi2e^''^PqrPAp + W^^e^'''^Pp,rPAp (AlO) 

and 

j^{cc)^Hu{cc) ^ 2£ief{vQRQ + VtRt + VttRtT cos 20 + VzzRzz + {VxzRxz - VOxB^x) COS 4) 

- VOzRqz + hiVyzRyz COS (j) - VOyRoy COS (f) - VxyRxy)}, (All) 

where 

Ryz COS 4> = i {W^/ - W'f) , (A12a) 

Roy cos (t) = i - ^a) ■ (Al2b) 

Note that the response functions R'y^ and R'^y are related to Wg-io terms which vanish in the plane- 
wave limit. It follows from the expressions (A7) and (All) that the cross sections asymmetry, which 
is measured at azimuthal angles = 7r/2 and (j) = — 7r/2, vanishes in the absence of the FSI. 
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